Quantum uncertainty relation saturated by the eigenstates of the harmonic oscillator 
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We re-derive the Schrodinger-Robertson uncertainty principle for the position and momentum of 
a quantum particle. Our derivation does not directly employ commutation relations, but works by 
reduction to an eigenvalue problem related to the harmonic oscillator, which can then be further 
exploited to find a larger class of constrained uncertainty relations. We derive an uncertainty relation 
under the constraint of a fixed degree of Gaussianity and prove that, remarkably, it is saturated 
by all eigenstates of the harmonic oscillator. This goes beyond the common knowledge that the 
(Gaussian) ground state of the harmonic oscillator saturates the uncertainty relation. 



The Heisenberg uncertainty principle [l[ captures the 
difference between classical and quantum states, and sets 
a limit on the precision of incompatible quantum mea- 
surements. It has been introduced in the early days 
of quantum mechanics, but its form has evolved with 
the understanding and formulation of quantum physics 
throughout the years. The first rigorous mathematical 
proof of Heisenberg's uncertainty relation for the canoni- 
cal operators of position x and momentum p ([x,p] = ih) 



((x 2 ) - (x) 2 )((p 2 ) (p) 2 ) > h 2 /4, 



(1) 



is due to Kennard [2j and Weyl [3| , but only pure states 
were considered there. The full proof was later derived 
following different methodologies 0-0] (see Q for more 
details), while the properties of the states saturating this 
inequality were also progressively unveiled. 

The original uncertainty relation ([T]) only concerned 
the operators x and p, but it was generalized to any pair 
of Hermitian operators by Schrodinger Q and Robertson 
H, in the case of pure states. In the same works, the 
anticommutator of x and p was also included in Eq. ([1} , 
yielding a stronger uncertainty relation 



((x 2 ) - (x) 2 )((p 2 ) (pf) 
-h(xp + px)-2(x) (p)f>h 2 /A, 



(2) 



that bears their name. The first proof of the Schrodinger- 
Robertson (SR) uncertainty relation for position and mo- 
mentum, Eq. ([2]), in the general case including mixed 
states is probably due to Moyal @, and, for any pair 
of not-necessarily Hermitian operators, to Dodonov, 
Kurmyshev and Man'ko [13] • In this latter work, the 
states of minimum uncertainty or minimizing states 
(MSs) for the SR inequality were identified as the pure 
states with a Gaussian wavefunction. Such Gaussian 
states are ubiquitous in physics as they play a major role 
for example in quantum optics (e.g., coherent states of 
the light field [ll|), atomic ph ysic s (e.g., collective exci- 
tations of an atomic ensemble jT2(]), optomechanics (e.g., 
nanomechanical oscillators fl3|). supraconductivity (e.g. 
superconducting LC circuits [lj]), etc. 

In this Rapid Communication, we revisit the status of 
Gaussian states in the context of uncertainty relations by 



exhibiting a connection with the harmonic oscillator and 
showing that remarkably, all its eigenstates - not just its 
ground state - appear as minimum uncertainty states. 
We first re-derive the SR inequality and corresponding 
MSs by using a variational method and standard alge- 
braic tools. This derivation reveals the direct link be- 
tween the quadratic order of Eqs. ([1]) and ([2]) in x and 
p and the fact that we deal with the quadratic Hamilto- 
nian of a harmonic oscillator. Then, we move on to find 
bounded uncertainty relations (To| . which give stronger 
bounds than Eq. ([2]) for states on which some a priori 
information is known, such as their purity Q or entropy 
|16|. Specifically, we derive a Gaussianity-bounded uncer- 
tainty relation, depending on the degree of Gaussianity 
of the state as measured by a parameter g that we intro- 
duce. We identify its corresponding set of MSs and find 
among them all the eigenstates of the harmonic oscilla- 
tor. This yields a fundamental new set of non-Gaussian 
minimum-uncertainty states, going beyond the common 
knowledge on the Heisenberg principle. 

Although the uncertainty relations, being at the root 
of quantum mechanics, have been investigated in various 
situations, such as multi-dimensional (rn . f]~sj ] or mixed 
states 0, EE HH, our results imply that there is more 
to gain by analyzing them under the perspective of the 
Gaussian character of a state. Non-Gaussian states of 
light can now be handled in the lab poM^ and have been 
proven essential in the field of continuous- variable quan- 
tum information [26l - [30l ] . but they remain hard to clas- 
sify. Identifying states of minimum uncertainty among 
them may lead to a better understanding of the struc- 
ture of the state space in infinite dimension and, since 
the Heisenberg principle is at the heart of the limitations 
on measurement precision |3l| . to the possible develop- 
ment of novel concepts in quantum metrology. 

Unconstrained SR relation. — We introduce our method 
as a way to find the MSs of the SR uncertainty relation. 
Consider a quantum state characterized by its density 
operator p. Its covariance matrix is defined as 
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| (xp + px) - (x) (p) 



\ {xp + px) - (x) (p) (p 2 ) - (p)' 



(3) 

where (•) = Tr(/5 •) stand for quantum expectation values 
in state p. Hereafter, we define the uncertainty of the 



2 



state p as the dimensionless variable 

a = 2(detj)^ 2 /h. (4) 

which is simply the square root of the left-hand side of 
Eq. ([2]) divided by H/2. We write it with the determi- 
nant of 7 to emphasize that it remains invariant under 
any linear canonical transformation in x and p, that is, 
under any operation modeled by a Hamiltonian that is 
quadratic in x and p (see Appendix A). This brings the 
strong simplification that it is sufficient to confine our 
search for MSs among states that satisfy the constraints 

Tr (px) = Tr (pp) = (5) 
Tr(p(xp + px)) = Tr (p (x 2 - p 2 )) =0 (6) 

that is, whose mean values vanish and 7 is proportional 
to the identity. Under these conditions, Eq. (Q} can be 
expressed as 

a = Tr05(l + 2n)), (7) 

where n — (x 2 +p 2 — l)/2 is the number operator for the 
harmonic oscillator (from now on, we assume h = 1 and 
take the particle's mass and angular frequency equal to 
one). 

Now, let us proceed with the minimization of a under 
constraints ([3]) and © by using the Lagrange multipliers 
method. For any density operator p, there is an eigenba- 
sis {|<£ n )} such that p = |*„) (* n | with < c„ < 1 
and c n — 1 ■ It is more convenient to define the unnor- 
malized vectors \tpn) — y/cn\^n) and rewrite the state 
as p = ^2 n \ip n ) (ip n \i while imposing the additional con- 
straint 

Tr(/5) = 1. (g) 

Then, choosing an orthonormal basis {\i)} to decompose 
the vectors \ip n ) — ^2 V'n l*)> we can re-express the uncer- 
tainty J7|) and constraints (JSJ) , ©, and flSJ) as functions 
of the Vn' s - We define the Lagrange multipliers X' k and 
consider the "uncertainty" functional 

a = a + A'jTr (p) + X' 2 Ti (px) + X' 3 Tr (pp) 
+ \' 4 Tr (p (xp + px)) + A' 5 Tr (p (x 2 - p 2 )) , (9) 

which implicitely depends on the complex amplitudes 
■0^'s. Extremizing a yields conditions on these ampli- 
tudes (see Appendix A) , which read as conditions on the 
unormalized eigenvectors \ip„) defining the minimizing 
state p, namely 

[ft + 1/2 + Ai + X 2 x + X 3 p + X 4 (xp + px) 

+\ 5 (x 2 -p 2 )]\Tjj n ) -0 (fO) 

where Afe = 2A' fc . Introducing the Hermitian operator 
H = n+l-+X 2 x+X 3 p+X 4 (xp + px)+X 5 (x 2 - p 2 ) , (11) 



we can rewrite Eq. (flQ)) as H \ip n ) = — Ai |V>n)i Vn, lead- 
ing to the necessary condition that the eigenvectors \ip n ) 
defining the MSs must be degenerate eigenvectors of H 
corresponding all to the same eigenvalue. 

Thus, one should diagonalize H in order to proceed 
with the identification of the MSs. As explained in Ap- 
pendix A, there exists a linear canonical transformation 
in x and p that transforms H onto the Hamiltonian of 
the harmonic oscillator, Ho = h + 1/2. Obviously, this 
means that H has the same eigenvalues as Hq and that its 
eigenvectors are U \n), where |n) are the number states 
(eigenstates of Hq) and U is the Gaussian unitary cor- 
responding to this canonical transformation. Since Hq 
does not possess any degeneracy in its spectrum, the 
same holds for H and therefore the only possibility is 
that the MSs is a pure state of the type p = U \n) (n\ U*. 
Among these states, we must keep those satisfying con- 
straints ([S])-©, which are the number states \n), so that 
the state that minimizes the uncertainty ([7]) is obviously 
the ground state |0). Of course, by plugging p = |0) (0| 
into Eq. (O, we recover the lower bound of the SR rela- 
tion, h 2 /A. By acting with linear canonical transforma- 
tions on |0), we obtain all Gaussian pure states, which is 
the well-known set of MSs for the SR uncertainty relation 
(see Appendix A). This was a long detour to re-derive 
Eq. © , but this connection with the harmonic oscillator 
turns out to be crucial in what follows. 

Degree of Gaussianity. — Our method works by reduc- 
tion to a constrained optimization problem (even for solv- 
ing the unconstrained SR inequality), so it can be sim- 
ply adapted to find the MSs with an extra constraint on 
Gaussianity. Several measures of non-Gaussianity have 
been used in the literature [H, [l(| HH-S], but we in- 
stead suggest using a parameter g capturing the degree 
of Gaussianity, inspired from our former work on non- 
Gaussian states with positive Wigner function (3f| [36| . 
Denoting as pa the Gaussian state that has the same co- 
variance matrix 7 (and same mean values (x) and (p)) as 
state p, we define the Gaussianity of p as 

g = Tr ( PPG ) /Tr (p 2 G ) . (12) 

It is more appropriate for our purposes and also has mer- 
its on its own, see Appendix B. Its main properties are as 
follows (see Appendix B for the proofs): (i) g is invariant 
under linear canonical transformations in x and p; (ii) g 
is a bounded quantity, that is, 2/e < g < 2, and g = 1 
for Gaussian states (but the converse is not true); (in) g 
provides a necessary criterion for the strict positivity of 
the Wigner function of a state. 

Let us briefly address possible experimental means for 
a direct estimation of g without going through a full state 
tomography procedure of p (see Appendix B). The trace 
overlap Tr (ppc) can be estimated using the eight-port 
homodyne detection scheme that is usually employed to 
estimate the Q-function of a state in quantum optics by 
simultaneously measuring the x and p quadratures of the 
two output modes of a balanced beam splitter [ll[. If 
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we inject the state p together with pa (instead of the 
vacuum state |0)) in this beam splitter preceding the 
quadrature measurements, then the value of the modified 
Q-function at the origin would read Q (0, 0) = Tr (ppc)- 
This method for measuring g would, however, require 
first performing homodyne measurements on p in order 
to obtain its covariance matrix 7 and prepare the Gaus- 
sian state pa- 

Gaussianity-bounded SR relation. — Our optimization 
method provides a necessary condition on the extremal 
solutions since it relies on the Lagrange multipliers 
method, so that concluding on a solution may become 
complicated if the eigenvectors of H cannot be identified 
analytically. Fortunately, finding an uncertainty relation 
under a fixed-g constraint leads to an eigenvalue problem 
that is analytically solvable. As before, we can confine 
our search on states satisfying constraints ([5]) and ([6|) 
since neither a nor g change under linear canonical trans- 
formations. For these states, the corresponding Gaussian 
state pa can be expressed simply as 



PG 



IN 



(13) 



where e 13 — ^4 and N 

a+l 



(a + l)/2. In addition, 

Tr (p G ) =l/a. Instead of minimizing the uncertainty a 
for a fixed Gaussianity g, it is easier to fix a and search for 
states of extremal g, or extremal overlap o = Tr (pe~^ n ) . 
The variational procedure for deriving the MSs is com- 
pletely analogous, but we extremize o using a constraint 
on the uncertainty in addition to Eqs. dS])-© and ©. 
Thus, the "overlap" functional is 

5 = o + X[Tr (p) + X' 2 Tr (px) + \' 3 Tr (pp) 
+ A 4 Tr (p (xp + px)) + X' 5 Tt (p (x 2 - p 2 )) 
+ A^Tr( / 5(2n+l)) (14) 

depending on six Lagrange multipliers A^. The extrem- 
ization conditions on the unnormalized eigenvectors \ipn) 
of the solution state can be written as 



-/3n 



+ 1/2+Ai + X 2 x + A 3 p + A 4 (xp + px) 



+ \ 5 (x 2 ~p 2 )+\ 6 n]\iP n )=0. (15) 
where A& = 2A' fc . By defining the Hermitian operator 
Hi =e~ 0fl + \ 2 x + \ 3 p + A 4 (xp + px) 



+ A 5 (x — p ) + A 6 n 



(16) 



we conclude that \ip n ) should be degenerate eigenvectors 
of Hi. It can be shown that, without loss of generality, 
we can restrict ourselves to states \ip n ) = ip n \n), with ip n 
being complex amplitudes (see Appendix B). For these 
states, the constraints (JS])-© are satisfied, and Hi is 
replaced by 



Ho 



- f3fl + A. 



6 n. 



(17) 



are possible if Xq < 0. As a result, we look for mixtures 
of two number states p — \ip n \ 2 \n) (n\ + \ip m \ 2 \ fn) (m\ 
satisfying the normalization constraint | ip n \ +\ip m \ =1 
and uncertainty constraint 

|Vg 2 (2n + 1) + |VAn| 2 (2m + 1) - a, (18) 
that achieve the minimum or maximum 

9=Wn\ ~ n+1 +\lpm\ - ■ ^ m+ i ■ (19) 



(a + l)' 



(a + l) 



By supervision, one can see that the minimum g (cor- 
responding to positive eigenvalues of H 2 ) is achieved by 
mixtures of two successive number states 



Pmin = r \n) (n| +(l-r)\n+ 1) (n + 1| . 



(20) 



with the parameters n and r € [0, 1[ depending on a. 
The number states \n) are naturally included in the set 
for r = 0. The maximum g (corresponding to negative 
eigenvalues of H 2 ) is achieved by mixtures 



Pmax = r |0) (0| + (1 - r) |n) (n\ 



(21) 



in the limit n -+ 00, r — > 1, while a = r + (1 — r)(2n + 1) 
is kept finite. 

In Figure [1] we plot as a dashed line the two extremal 
values of g for a fixed a as realized by the states of 
Eqs. ([2"0 |) -([2"T ]) . The MSs (i.e., the states minimizing the 
uncertainty a for fixed degree of Gaussianity g) corre- 
spond only to some part of this line, which we show as 
a solid line. For g > 1, the situation is simple and all 
states p m ax are MSs. For this region, the minimum on 
the uncertainty as a function of g can be easily derived 



a > 
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2-5' 



if g > 1, 



(22) 



employing Eq.(21) (see Appendix D). In contrast, for - < 
g < 1, the minimum a for fixed g displays discontinuities. 
In the interval 



(n + l) n+1 (3 + 2n) n n (1 

2^ " <9< 



2n) 



(2- 



(1 



(23) 



a is minimized by the states /5 m i n with the value of n 
satisfying Eq. (f2"3"]l . In particular, we see that all number 
states \n) are included in this set for specific values of g 
corresponding to the upper bound of Eq. (|23[) for different 
n's. Moreover, in Fig. [TJb), we see that most of the 
MSs (solid line) consist of these states \n) and their close 
neighborhood. Thus, as advertised, we conclude that all 
eigenstates of the harmonic oscillator are extremal in the 
sense that they exhibit the lowest allowed uncertainty 
given their non- Gaussian character. 

For other values of g, once the value of n is identified 
from Eq. (|23[) . one has to solve the polynomial equation 



The eigenvectors of H 2 are the number states \n), but, 
unlike for the harmonic oscillator, double degeneracies 
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FIG. 1: (a) Extremal values of the degree of Gaussianity g 
for a fixed uncertainty a shown as a dashed (blue) line. It is 
achieved by states p m i n for g < 1, and p ma x for g > 1. The 
line connecting subsequent number states \n) and |n + 1) is 
realized by mixtures of them. The Gaussianity-bounded SR 
relation corresponds to the part of this extremal line shown 
as a solid (black) line. The uncertainty a must have a value 
larger or equal to the solid line for a given g. (b) Magnified 
view of figure (a), where the discontinuity of the uncertainty 
relation becomes evident. 



Pmin and /5 max , which are phase invariant, the set of MSs 
comprises all states with a covariance matrix propor- 
tional to the identity which can be transformed onto them 
by phase averaging. Furthermore, all states connected to 
p m j n and /5 max by linear canonical transformations are 
MSs as well, since the uncertainty a and Gaussianity g 
are invariants of the group. 

Conclusions. — We have exhibited a variational method 
to derive the Schrodinger- Robertson uncertainty relation 
by casting it as an eigenvalue problem related to the har- 
monic oscillator. It follows an "inverse path" to the com- 
mon procedure where the lower bound on the uncertainty 
is derived based on commutators, and the MSs are then 
identified. Such an inverse procedure was put forward 
by Dodonov and Man'ko for the derivation of purity- 
bounded uncertainty relations 0, [H| , but it appears that 
our method is more generally applicable because it is 
based on the amplitudes of the eigenstates of the MSs 
instead of its density matrix elements, see Appendix D. 
It is especially useful when constraints are included that 
account for some knowledge on the state. 

In particular, we have found a new uncertainty rela- 
tion that is bounded by the degree of Gaussianity g of 
the state. The state with the overall lowest uncertainty 
a is of course the Gaussian ground state |0) of the har- 
monic oscillator (g = 1), but we have thus proven that 
the non-Gaussian states with the lowest uncertainty a 
for a fixed g < 1 include as well mixtures of subsequent 
number states |n) of the harmonic oscillator. Among 
these MSs, the number states play a prominant role 
as they are the only phase-invariant pure states. We 
have proven that the number states are also extremal in 
this uncertainty-related sense, thereby extending to all 
(non-Gaussian) eigenstates of the harmonic oscillator the 
celebrated minimum-uncertainty property of its (Gaus- 
sian) ground state. Given the considerable attention that 
non-Gaussian states are attracting in continuous- variable 
quantum information theory, see e.g. [37] , unveiling this 
cxtremality property of harmonic oscillator states may 
contribute to further fundamental progress in the field, 
especially in relation with quantum metrology. 

AM gratefully acknowledges financial support from the 
F.R.S-FNRS. This work was also carried out with the fi- 
nancial support of the F.R.S-FNRS via project HIPER- 
COM and the support of the Belgian Federal program 
IUAP via project Photonics@be. 



for a in order to find the dependence of the uncertainty 
a on the Gaussianity g. For 3/4 < g < 1, which covers 
most of the interesting region, the explicit expression (see 
Appendix C) is 

a> 2 + 2 ^~V-l , if 3/4 <«?<!. (25) 
9 



Appendix A: Full proof of the SR uncertainty 
relation 

Consider an arbitrary quantum state defined by the 
density operator p. Using standard notations in quantum 
optics, the covariance matrix 7 of this state is 



In the Appendix B we show that, in addition to the states 



H^^iiifi-diU^-d^Yp) (Al) 
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where r = (x,p) T is the vector of position and momen- 
tum observables, d =Tr(rp) is the displacement vector, 
and {•, ■} stands for the anticommutator. We start from 
the definition of the dimensionless uncertainty 

(det 7 ) 1 / 2 
a = —hj^- (A2) 

The determinant of 7, and consequently the uncertainty 
a, remains invariant under the action of the linear canon- 
ical group, the semidirect product of the special linear 
group Sp(2, R) with thetranslation group T (2). In quan- 
tum optics, these correspond to Gaussian operations, 
combining displacements and symplectic transformations 
[TtJ ■ Therefore, without loss of generality, we may confine 
our search of MSs among states with a covariance matrix 
7 proportional to the identity (the so-called Williamson 
normal form) and d = 0. In other words, we may search 
for states p that minimize the uncertainty a while satis- 
fying the constraints 

Tr (px) = Tr (pp) = (A3) 
Tr (p(xp + px)) = Tr(p(x 2 -p 2 )) =0 (A4) 

Under these conditions, the uncertainty a can be ex- 
pressed as 

a = Tr (p (2n + 1)) , (A5) 

where n = a)a the number operator. 

Now, let us provide more details on the optimiza- 
tion of a under the constraints of Eqs. (IA3[) - (IA4|) . It 
is known that for every density matrix p, a unique 
eigenbasis {|\&n)} exists such that p = J2c n \^ n ) (^ n \ 
with < c n < 1 and Y] Cn = 1, It is more conve- 
nient for our purposes to define the unnormalized vectors 
\tpn) = \^\^n), and rewrite p as 

P = X>»)W»I (A6) 



while imposing the additional constraint 

Tr(p) = l. (A7) 

Then, we choose an orthonormal basis {\i)} to decompose 
the vectors \ip n ) — I*) an d accordingly re-express 

the uncertainty (|A5|) as 

a=J2i>n*Pi(i\(2n + V\j) (A8) 



and the constraints (|A3|) . (|A4p . and (fATf as 



n,i,j 


= 


(A9) 


\ y 1 i* / 7 / ■ 1 ~ i -\ 


= 


(A10) 


EC^n («|(ip + pi) |j) 


= 


(All) 


E« (*I(£ 2 -p 2 ) li) 


= 


(A12) 




= 1. 


(A13) 
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We define the Lagrange multipliers A' fc and consider the 
"uncertainty" functional 

5 = a + A^Tr (p) + A^Tr (/5.x) + A 3 Tr (pp) 
+ A^Tr (p (ip + pi)) + A^Tr (p (i 2 - p 2 )) , (A14) 

which implicitely depends on the complex amplitudes 
■0^'s via the relations Eqs. (|A9[ )- ([A13p . Extremizing a 
yields conditions on these amplitudes which read as con- 
ditions on the unnormalized eigenvectors \ip n ) defining 
the minimizing state p, namely 

[n + 1/2 + Ai + X 2 x + A 3 p + \ 4 (xp + px) 

+\ 5 (x 2 ~p 2 )}\^j n ) = (A15) 

where A^ = 2AJ,. Introducing the Hermitian operator 

H = h + - + X 2 x + A 3 p + A 4 (xp + px) + A 5 (x 2 - p 2 ) , 

(A16) 

we can rewrite Eq. (|A15|) as H \ijj n ) = — Ai \ip n ), Vra, lead- 
ing to the necessary condition that the eigenvectors \tp n ) 
defining the MSs must be degenerate eigenvectors of H 
corresponding all to the same eigenvalue. 

One should identify the eigenvectors and eigenvalues 
of H. In the language of quantum optics, this is an 
easy task because one can always apply a squeezing 
and displacement operation that transforms H onto the 
Hamiltonian of the harmonic oscillator, Hq = h + 1/2. 
More precisely, there exists a linear canonical transfor- 
mation U — exp (70^ — 7*a) exp (/3a^ 2 — f3*a 2 ) such that 
H = UHqW, up to a real constant. Obviously, this 
means that H has the same eigenvalues as Hq, up to this 
constant, and that its eigenvectors are the accordingly 
transformed number states, U \n), remembering that the 
number states \n) are the eigenstates of Hq. Among these 
states, we must keep those satisfying constraints (|A3[1 - 
(|A4[) . which are the number states \n), so that the state 
that minimizes the uncertainty (|A5I) is obviously the vac- 
uum state |0). 

Clearly the vacuum state |0) is the minimizing state 
just for the specific class of states that we have chosen for 
our proof i.e. the states with covariance matrix propor- 
tional to identity. The whole set of MSs can be obtained 
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by acting the linear canonical transformations on the |0) 
since our proof is invariant under the action of this group. 
In other words, the "orbit" of states that are connected 
to |0) by linear canonical transformations (that is, all 
Gaussian pure states) coincide with the well-known set 
of MSs for the SR uncertainty relation. 

Appendix B: Degree of Gaussianity 

In our work, we have introduced the quantity 

g = Tr(pp G )/Tr(p 2 G ) (Bl) 

in order to characterize the degree of Gaussianity of a 
state p. In Eq. (|B1|) , pa is the reference Gaussian state of 
p in the sense that it possess the same covariance matrix 
7 as p. Here, we exhibit the proofs of the mathematical 
properties of g that we list in the main text and elaborate 
on the physical intuition behind its definition. 

(i) The Gaussianity is invariant under linear canon- 
ical transformations, i.e., symplectic transformations 
Sp(2,R) and translations T(2). 

Proof. A Gaussian operation U G acting on a state p 
can be always translated as a sequence of displacement, 
rotation, and squeezing of the Wigner function of p in the 
phase-space [38|]. The Wigner function of the reference 
Gaussian state pa experiences the same deformation in 
phase space as that of p, and there ore one may con- 
clude that pq experiences the same Gaussian operation 
U G as p. Then, using the invariance of trace under cyclic 
permutations, we have 

Tr (lJ G pU G U GPG U G ) = Tr (pp G ) 

Tr (u G p G Ul U G p G U G ) = Tr (p 2 G ) , (B2) 

so we conclude that g is invariant under Gaussian oper- 
ations. 

(it) The Gaussianity is a bounded quantity, that is, 
2/e < g < 2. For Gaussian states g = 1, while the 
converse is not necessarily true. 

Proof. Having identified in this work the extremum 
values for g as a function of the uncertainty a, it is 
straightforward to derive the lower and upper bounds 
of the interval [2/e, 2], which correspond to the limit 
a — > oo. The state which realizes the lower bound (up 
to Gaussian transformations) is the number state \n) as 
n — > oo, while the upper bound is saturated by the state 

p max = (l-r)\0){0\ + r\n)(n\ (B3) 

with r — > and n — > oo. It is also straightforward to 
prove that if a state is Gaussian, then g = 1. On the 
other hand, g is not constructed on the basis of a mathe- 
matical distance, unlike the measures of non- Gaussianity 
of refs. |32h34{ . As a consequence, the inverse statement 
is not necessarily true and non-Gaussian states might ex- 
ist possessing a degree of Gaussianity equal to 1. 
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FIG. 2: Fhe red dotted line represents Eqs. (|B4|) . Fhe black 
solid line represents the bounds on Gaussianity g derived in 
the main text. No quantum state with completely positive 
Wigner function can exist in the area outside the two dotted 
lines. No quantum state can exist outside the two solid lines. 

(Hi) The Gaussianity provides necessary criterion for 
strict positivity of the Wigner function of a state. 

Proof. In a previous work (3(| , we have derived bounds 
on the trace overlap Tr {pp G ) for states with strictly pos- 
itive Wigner function. These bounds have been derived 
partially analytically and can be easily translated into 
bounds on the Gaussianity g. Based on the formulas de- 
rived in [3(1 ] , we can prove that a state with strictly posi- 
tive Wigner function and of uncertainty a is restricted to 
have a Gaussianity g that lies in the interval [g-mm , Smax] , 
where 

5min = 0.0095a + 0.62 + 0.711/a - 0.333/a 2 



5max=J— — . (B4) 
V 1 + 5* 

It should be noted that the bounds on the Gaussian- 
ity g given by our Gaussianity-bounded uncertainty rela- 
tion provide a necessary criterion for a quantum state to 
be physical. The classical bivariate distributions which 
achieve the minimum degree of Gaussianity among all 
continuous positive-definite distributions (classical and 
quantum quasi-distributions) of the same covariance ma- 
trix have been identified in a previous work (3o| . Inter- 
estingly, one can always find a classical distribution that 
possess a smaller g than what is allowed for all quan- 
tum states. This indicates that the positivity of the den- 
sity matrix that is imposed to derive the Gaussianity- 
bounded uncertainty relation is more restrictive than the 
positivity and continuity of a distribution that is imposed 
on the proof of [35]]. In FigJ2]we represent the Eqs. (IB4) 
by dotted red lines, while the ultimate bounds on g, as 
these derived in the main text, by solid black lines. For 
a state outside the dotted lines we know that its Wigner 
fuction has necessarily negative parts. Outside the black 
lines no quantum states can exist. 
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To obtain more intuition about the quantity g, it is 
instructive to explicitly express it in terms of the mo- 
ments of the Wigner quasi-probability distribution. Let 
us restrict ourselves to a state p with a covariance ma- 
trix 7=5 diag(a, a), since by linear canonical transfor- 
mation all states can be reduced to this form and we 
have proven that g remains invariant along this transfor- 
mation. In the Wigner representation, the correspond- 
ing reference Gaussian state pa of such a state is phase- 
invariant; hence, it can be expressed as 



W G (r) = — e 
ira 



-r 2 /c 



\Jx 2 + 



p 2 . 



(B5) 



In contrast, the state p itself may possess an angular- 
dependent Wigner function W (r, ip) . The trace overlap 
between p and pa takes the following form in the Wigner 
representation 



Tr (pp G ) = 2ir / / W (r, p) W G (r) rdrdip 



(B6) 



while we have g — aTr (ppc)- 

Now let us average the phase of the Wigner function 
W (r, tp) in order to construct a new state p s with phase- 
invariant Wigner function W s (r), 



W s {r) 



2vr 



W (r, ip) dip. 



(B7) 



Employing Eq. (IB7[) . the trace overlap Eq. (IB6|) can be 
re-written as Tr (ppc) = 4ir 2 J W s (r) Wq (r) rdr, and by 
expanding W G (t) in Taylor series we arrive to 



47T 



oo 

^ [ n\ (a) n 



(B8) 



For a Gaussian bivariatc distribution independent of the 
phase as in Eq. (|B5I) , a simple expression exists for the ra- 
dial moments (?" 2n+1 ) G = {oi) n T (n + 1) /2n. Since by 
definition (r 3 ) = (r 3 ) G , we conclude that the Gaussian- 
ity g accounts for the difference of the odd (> 5) radial 
moments of the Wigner function of the phase-averaged 
state p s as compared to those of the reference Gaussian 
state. 

We can use now the phase averaged state p s intro- 
duced in Eq. (|B7p to prove that in our search for quan- 
tum states which extremize g while possessing a covari- 
ance matrix of the form 7 = ■§ diag(a, a), we are allowed 
to restrict to mixtures of number states. By substitution 
of Eq. jB7]) into Eq. ([B6]) . it is evident that 



Tr (pp G ) = Tr {p s p G ) 



(B9) 



In addition the reference Gaussian state (and therefore 
the covariance matrix) of p s is the same as for p, since 
phase averaging cannot affect the phase-independent 
Wigner function Eq. (|B5j) . From this, one may conclude 
that g is the same for p and p s and therefore, invariant 
under the phase randomization procedure for states with 



local oscillator 




FIG. 3: The modified eight-port homodyne scheme which can 
be used to estimate Gaussianity g. In the usual scheme which 
used to estimate the Q-function of a state p, the vacuum state 
|0) is fed in the down-right beam splitter instead of pa- 



covariance matrix of the form 7 = -| diag(a, a). There- 
fore, we may with no loss of generality confine our search 
to states with phase-independent Wigner function, which 
can be always expressed as mixtures of number states, 
i-e-, p = J2°n \n) (n\. 

Finally, as explained in the main text, there are pos- 
sible experimental means for a direct estimation of the 
Gaussianity g of a state p, without going through a full 
state tomography procedure. The trace overlap between 
p and pa can be estimated without the need for a full 
state tomography, usin g th e modified eight-port homo- 



dyne detection scheme 



|39f shown in Fig[3] However, 
this would first require performing homodyne measure- 
ments on p for the identification of its covariance and 
the reconstruction of the reference Gaussian state pg- 
As a consequence, the whole procedure requires multiple 
copies of p just as does a full state tomography procedure, 
but on the other hand it avoids the instabilities which ap- 
pear in the reconstruction of the Wigner function from 
experimental data. Other scenarios for measuring g with- 
out constructing pg from p may exist, and it is a subject 
that deserves further investigation. 



Appendix C: Analytic expression for the 
Gaussianity-bounded uncertainty relation 

From the general solution for the Gaussianity-bounded 
MSs as derived in the main text, we notice that for the 
interval | < g < 1 the uncertainty a is saturated by 
mixtures of |0) and |1). The dependence of the minimum 
value of the uncertainty on Gaussianity degree can be 
easily derived by solving the equation 



4a- 



(<*-!)" 
(a + 1) 2+ " 



(1 + n) =9 



(CI) 



for a when n — 0. One arrives to the following inequality 



a > 



2 + 2VT^ff-.g 



if i<5 <l. 
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For g > 1, the MSs have been proven in the main text 
to comprise the following set, 



p=(l-r)|0) (0|+r» (n\ 



(C2) 



in the limit n — > oo, r — ¥ 0, and while a = (1 — r) 
r(2n + 1) is kept finite. Employing the identity 



(1-r) 



2a 



(a + 1) (a+1) 



r 2 "^". 1 ?!* (03) 



71+1 



one can derive the following dependece of the lower limit 
of uncertainty as a function of Gaussianity, 



a > 



fJ 



9 



if g > 1. 



Thus, ignoring the remaining tiny interval - < g < § 



(note 



0.014), we can summarize the Gaussianity- 



bounded uncertainty relation as 

a > 



„ - if | < 5 < 1 



if 5>!- 



Gaussianity-bounded uncertainty relation, is applicable 
not only in the case where all constraints are linear in 
the density matrix elements of p, i.e., of the form Tr (Bp), 
but also in the case where we have non-linear constraints 
such as the purity Tr (p 2 ) or the von Neumann entropy 
— Tr(/51n/5). It is not difficult to show that in the pres- 
ence of non-linear constraints, the necessary condition 
on the existence of degeneracies in the spectrum of the 
Hermitian operator constructed from the constraints is 
lifted. The necessary condition is, in this case, that ev- 
ery eigenvector \ij) n ) of the density matrix of the solution 
p = IV'n) (V'nl should be an eigenvector with positive 
eigenvalue of a Hermitian operator derived in a similar 
way as in the main paper 



(F-p)|V„} = 0. 



(Dl) 



Moreover, as it is dictated by the above condition 
Eq. (|D1[) . the eigenvalues are the mixing amplitudes of 
the eigenvectors 



H\lp n ) = {lp n \ 4>n) \lpn) = C n \lpn) 



(D2) 



Appendix D: Generalization to non-linear 
constraints 

The method we have developed here in order to 
re-derive the SR uncertainty relation and derive the 



In a further work [40j , we explicitly show via an example 
how the method develops in this obviously more com- 
plicated case and we draw parallels with the method 
of derivation of purity -bounded uncertainty relations 
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